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Abstract: The article deals with electrodynamics in the presence of anisotropic 
materials having scalar wave impedance. Maxwell's equations written for differ- 
ential forms over a 3-manifold are analysed. The system is extended to a Dirac 
type first order elliptic system on the Grassmannian bundle over the manifold. 
The second part of the article deals with the dynamical inverse boundary value 
problem of determining the electromagnetic material parameters from bound- 
ary measurements. By using the boundary control method, it is proved that the 
dynamical boundary data determines the electromagnetic travel time metric as 
well as the scalar wave impedance on the manifold. This invariant result leads 
also to a complete characterization of the non-uniqueness of the corresponding 
inverse problem in bounded domains of M"^ . 



Introduction 

Classically, the laws of electromagnetism expressed by Maxwell's equations are 
written for vector fields representing the electric and magnetic fields. However, 
it is possible to rephrase these equations in terms of differential forms. It turns 
out that this alternative formulation has several advantages both from the the- 
oretical and practical point of view. First, the formulation of electromagnetics 
with differential forms reflect the way in which the fields are actually observed. 
For instance, flux quantities are expressed as 2-forms while field quantities that 
correspond to forces are naturally written as 1-forms. This point of view has 
been adopted in modern physics at least when fields in free space are dealt with, 
see Furthermore, the formulation distinguishes the topological properties 
of the electromagnetic media from those that depend on geometry. It is under- 
stood that geometry is related to the properties of the material where the waves 
propagate. The distinction between non-geometric and geometric properties 
has consequences also to the numerical treatment of the equations by so called 
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Whitney forms. An extensive treatment of this topic can be found in jl^], |jT^. 
For the original reference concerning Whitney elements see ]58[ |. 

The present work is divided in two parts. In the first part, we pursue further the 
invariant formulation of Maxwell's equations to model the wave propagation in 
certain anisotropic materials. More precisely, we consider anisotropic materials 
with scalar wave impedance. Physically, scalar wave impedance is tantamount 
to a single propagation speed of waves with different polarization. The invari- 
ant approach leads us to formulate Maxwell's equations on 3-manifolds as a first 
order Dirac type system. From the operator theoretic point of view, this formu- 
lation is based on an elliptization procedure by extending Maxwell's equations 
to a Grassmannian bundle over the manifold. This is a generalization of the 
elliptization of Birman and Solomyak and Picard (sce[p], [^6| ). 

In the second part of the work, we consider the inverse boundary value problem 
for Maxwell's equations. In terms of physics, the goal is to determine material 
parameter tensors, electric permittivity e and magnetic permeability fi, in a 
bounded domain from field observations at the boundary of that domain. As 
it is already well established, for anisotropic inverse problems it is natural to 
consider the problem in two parts. First, we consider the invariant problem on 
a Riemannian manifold, where we recover the travel time metric and the wave 
impedance on the manifold. As a second step, we consider the consequences of 
the invariant result when the manifold is imbedded to M'^. 

Although inverse problems in electrodynamics have a great significance in physics 
and applications, results concerning the multidimensional inverse problems are 
relatively recent. One-dimensional results have existed starting from the 30'ies, 



see e.g. [Q, |50 . The first breakthrough in multidimensional inverse problems 
for electrodynamics was based on the use of complex geometrical optics | |52| , 
[ p^ , p3[ , iQ. In these papers, the inverse problem of recovering the scalar ma- 
terial parameters from complete fixed frequency boundary data was solved even 
in the non-selfadjoint case, i.e., in the presence of electric conductivity. These 
works were based on ideas previously developed in references |]5^ , , | to 



solve the scalar Calderon problem, that obtained its present formulation in Hi'? 

In the dynamical method to solve an isotropic inverse boundary problem 

based on ideas of integral geometry is developed in [Q. The method, however, 
is confined to the case of a geodesically simple manifolds and, at the moment, 
is limited to finding some combinations of material parameters, including elec- 
tric conductivity. An alternative method to tackle the inverse boundary value 
problem is the boundary control (BC) method, originated in [Q. Later, this 
method was developed for the Laplacian on Riemannian manifolds and for 
anisotropic self-adjoint - |2|] and certain non-selfadjoint inverse problems 
[32|. The first application of the BC method to electrodynamics was done in 
[^T @]' authors of these articles show that, when the material parameters 
e and ^ are real scalars or alternatively when e — fi, the boundary data deter- 
mines the wave speed in the vicinity of the boundary. These works employed 
the Hodge- Weyl decomposition in the domain of influence near the boundary. 
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The real obstruction for this technique is that, as time grows, the domain of 
influence can become non-smooth and the topology may be highly involved. For 
these reasons, our paper is based on different ideas. 

In this article, there are essentially two new leading ideas. First, we character- 
ize the subspaces controlled from the boundary by duality, thus avoiding the 
difficulties arising from the complicated topology of the domain of influence. 
The second idea is to develop a method of waves focusing at a single point of 
the manifold. This enables us to recover pointwise values of the waves on the 
manifold. The geometric techniques of the paper are presented in po| and the 
book ||. 

The main results of this paper can be summarized as follows. 

1. The knowledge of the complete dynamical boundary data over a suffi- 
ciently large finite period of time determines uniquely the compact man- 
ifold endowed with the electromagnetic travel time metric as well as the 
scalar wave impedance (Theorem 4.1). 

2. For the corresponding anisotropic inverse boundary value problem with 
scalar wave impedance for bounded domains in M.^ , the non- uniqueness is 
completely characterized by describing the class of possible transforma- 
tions between material tensors that are indistinguishable from the bound- 
ary (Theorem 11.1). 

To the best knowledge of the authors, no global uniqueness results for inverse 
problems for systems with anisotropic coefficients have been previously known. 
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to thank Helsinki University of Technology, Loughborough University and Uni- 
versity of Helsinki for their kind hospitality and financial support. Furthermore, 
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1 Maxwell's equations for forms 

In this chapter we derive an invariant form for Maxwell equations, consider 
initial boundary value problem for them and show how energy of fields can be 
found using boundary measurements. 

We start with Maxwell equations in domain C M'^ equipped with the standard 
Euclidean structure. Since our objective is to write Maxwell equations in an 
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invariant form, we generalize the setting in very beginning and instead of domain 
r2 consider manifolds. 

Let (M, (7o) be a connected, oriented Riemannian 3-manifold possibly with a 
boundary dM ^ 0. We assume that all objects in this paper are C°°-smooth. 
Consider Maxwell's equations on Af , 

curls = -Bt, (Maxwell-Faraday), (1) 
curliJ = Dt, (Maxwell-Ampere), (2) 

where E and H are the electric and magnetic fields, and B and D are the 
magnetic flux density and electric displacement, assumed for the time being to 
be smooth mappings M x R ^ TM . Here TM denotes the tangent bundle 
over M . The curl operator as well as divergence appearing later will be defined 
invariantly in formula ^ below. The sub-index t in the equations (|l|)-(H) 
denotes differentiation with respect to time. We denote the collection of these 
vector fields as T{M x R). At this point, we do not specify the initial and 
boundary values. To avoid non-physical static solutions, the above equations 
are augmented with the conditions 

divS ^ 0, divD = 0. (3) 

Furthermore, the fields E and Z?, and similarly the fields H and B are inter- 
related through the constitutive relations. In anisotropic and non-dispersive 
medium, the constitutive relations assume the simple form 

D = eE, B = ^iH, (4) 

where e, are smooth and strictly positive definite tensor fields of type (1, 1) on 
M . Our aim is to write the above equations using differential forms. 

Given the metric (70, we can associate in a canonical way a differential 1-form to 
correspond each vector field. Let us denote by a'^T*M the fc:th exterior power 
of the cotangent bundle. We define the mapping 

TM T*M, X ^ 

through the formula go{X,Y) — X^{Y). This mapping is one-to-one and it 
has the following well-known properties (See e.g. ||5^): For a scalar field u G 
C°°{M), (gradu)'' — du, where d is the exterior differential and for a vector 
field X G r(M), we have 

(curlX)^ = *o^^^^ divX = ~SoX\ (5) 

where *o denotes the Hodge-* operator with respect to the metric go, 

*o : a''T*M -> A^-''T*M, 
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and Sq denotes the codifFerential |^ 

So = (-1)'= *o d*o : n^M n^-^M. 

Here, fi'^M denotes the smooth sections M — > a'^T*M, i.e. differential /c— forms. 
Applying now the operator b on Maxwell's equations yields 

dE^ = -*oB\, dH^^*QD\, 

where we used the identity *o*o = id valid in 3-geometry0. The divergence 
equations (||) read 

5oD^ = 0, 5oB^ = 0. 

Consider now the constitutive relations (^. Starting with the equation D = eE, 
we pose the following question: Is it possible to find a metric such that the 
Hodge-* operator with respect to this metric, denoted by *£, would satisfy the 
identity 

Assume that such a metric exists. By writing out the above formula in 
given local coordinates (x^, a;^, x"^) and recalling the definition of the Hodge-* 
operator, the left side yields 

*o(e-B)^ = ^g^ejpqgo^.je'^^E^dxP A dx« 



where e is the totally antisymmetric permutation index and go = det((7o,jj)- 
Likewise, the right side reads 

*^E^ = ^^gi^ejpqgo^^kE'^dxP A 

so evidently the desired equality ensues if we set 

By taking determinants of both sides we find that 

Vg'e = Vgodet(e). 

Thus we see that the appropriate form for the metric tensor in the contravariant 
form is 

In the same fashion, we find a metric such that 



^Cf. with <5o = (—1)"*+"+! *g cl*Q for Riemannian n-manifolds 
^For Riemannian n-manifold, we have in general *o*o = (— l)'''""'^' 
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In general, the metrics and can be very different from each other. In this 
article, we consider a particular case. Indeed, assume that the material has a 
scalar wave impedance, i.e., the tensors e and fi satisfy 

2 

^ = a e, 

where the wave impedance, a = a{x), is a smooth function on M. Now we 
define two families of 1- and 2-forms on M as follows. We set 

Similarly, we define 

u^^aH\ i^^ = *oaD\ (7) 

Observe that the wave impedance scaling renders and r]^ to have the same 
physical dimensions, and the same holds for the 2-form. Now it is a straight- 
forward matter to check that the constitutive relations assume the form 

2 12-^1 

a 

We can make these equations even more symmetric by proper scaling of the 
metrics. Indeed, since a~^/x = ae, we have a new metric g that is defined as 



9 = 9ae 



We have, by direct substitution that 

This new metric will be called the travel time metric in the sequel. 
Assume that 

denotes the Hodge-* operator with respect to some metric g. If we perform a 
scaling of the metric as 

the corresponding Hodge operator is scaled as 

Therefore, if we denote by * the Hodge-* operator with respect to the travel 
time metric, we have 

* = a*e = -*u : A^T*M a'^T*M. 

a 

But this means simply that, in terms of the travel time metric, we have 

i^^ = *uj\ = *i/^ (9) 



6 



Consider now Maxwell's equations for these forms. After eliminating the v- 
forms using the constitutive equations (|^), Maxwell-Faraday and Maxwell- 
Ampere equations assume the form 



I Saio^ = Lul, Sa = (-1)'" * ad-* : Q'^M Vl'^'^M (10) 



and the divergence equations (H) read 



du^ 



0, 5a,UJ^ = 0. 



(11) 



In the sequel, we call equations (jT^) and ( pi] ) Maxwell's equations. 

It turns out to be useful to define auxiliary forms that vanish in the electromag- 
netic theory. Let us introduce the auxiliary forms oj*^ and via the formulas 



Furthermore, we define the corresponding i^-forms as 



(12) 



Since these auxiliary forms are all vanishing, we may modify the equations ( po|) 
to have 

dio^ — SaiiJ^ — —i^tJ dioo ~ <5qW^ = ~'^t- (1>^) 

Putting the obtained equations together in a matrix form, we arrive at the 
equation 

ujt+MLu = 0, (14) 



where 



and the operator A4 (without defining its domain at this point, i.e., defined as 
a differential expression) is given as 



M 



( 


-5a 
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d 





-5a 








d 





-5a 




V 





d 








(15) 



The equation ( p4| ) is called the complete Maxwell system. In the next section, 
we treat more systematically this operator. 

Remark 1. The operator M has the property 

= -diag(A°, Ai, A^, A^) = -A„, 
where the operator A^ acting on fc-forms is 



Al^d5a + 5ad^ A^g+Q{x,D), 
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with Ag denoting the Laplace-Beltrami operator on fc~forms with respect to 
the travel time metric and (5(x, D) being a first order perturbation. Hence, if 
Lu satisfies the equation (|lj), we have 

+ A„)lj = {dt - M){dt + M)lj = 0. 

In particular, we observe that the assumption that the impedance is scalar 
implies a unique propagation speed for the system. 

Remark 2. Denote by QM = ©l^gfi'^Af the Grassmannian algebra of differ- 
ential forms, where fi'^M are the differential fc-forms. Then the operator Ai in 
formula (|l5| ) can be also considered as a Dirac operator d — Sa '■ flM. 

Before leaving this section, let us briefly consider the energy integrals in terms 
of the differential forms. In terms of the vector fields, the energy of the electric 
field at a given moment t is obtained as the integral 

£{E)^ [ eE-EdV= [ go{E,D)dV = [ E^ A *oD^ 
Jm Jm Jm 

where dV is volume form of (Af, go). By plugging in the defined forms we arrive 
at 

£{E) = ( -Lo^ A*uj^. 
Jm oi 

In the same fashion, we find that the energy of the magnetic field reads 

£{B) = I -u? A ^u?. 
hi " 

These formulas serve as a motivation for our definition of the inner product in 
the following section. 

1.1 Maxwell operator 

In this section we establish a number of notational conventions and definitions 
concerning the differential forms used in this work. 

We define the L^-inner products for fc-forms in as 

Jm a 

Further, we denote by L'^{VL^M) the completion of with respect to the 

norm defined by the above inner products. We also define 

L2(M) = L^{n°M) X L^in^M) X L^{9?M) x L^{n^M). 

Similarly, we define Sobolev spaces H'*(M), s e M, 

W{M) ^ H'{n°M) X Wiyi^M) X H'in'^M) X H^'in^M), 
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H^(M) = H^{n°M) X H^ifl'^M) X H^ifi^M) x H^{n^M), 

where H'^{fl'^M) are Sobolev spaces of k- forms. At last, H^{n^M) is the 
closure in H''{Vt^M) of i.e. the subspace of of fc— forms which 

vanish near dM . 

The domain of the exterior derivative d in the L^-space of fc-forms is 

H{d,n''M) = {lo'' e L^(ifM) I dw'' e L^{^t+^M)] . 
Similarly, we set 

where 5a is the weak extension of the operator 5a ■ il'^M il''~^M. In the 
sequel, we shall drop the sub-index a from the codifferential. 

The codifferentiation 5 is adjoint to the exterior derivative in the sense that for 
C^-forms on M, 

To extend the adjoint formula for less regular forms, let us first fix some nota- 
tions. For uj^ e fl'^M, we define the tangential and normal boundary data at 
DM as 

a 

respectively, where i* : fl'^M Vl^dM is the pull-back of the natural imbedding 
i : dM M. Sometimes, we denote n = Hq, to indicate a particular choice a. 
With these notations, let us write 

JdM a 

We add here a small caveat that the above formula does not define an inner 
product as w*"' and 77*^+^ are differential forms of different order. For ut'^ £ Vl^M 
and g fi'^+^M, the Stokes formula for forms can be written as 

{duj\ r;'=+i)i. ™ (c.^ <57y^+i)i. = {tu\nrj^+^). (16) 

This formula allows the extension of the boundary trace operators t and n 
to H{d,n^M) and H{5,n^M), respectively. Indeed, ii uj^ e H^{n^M), then 
tw'= G H^/^{n''dM) and, by formula ®, we may extend 

t : H{d,n''M) H-^/'^{n^dM). 

In the same way, equation (|l6| ) gives us the natural extension 

n : H{5,n''+^M) = H-^/^{n^^''dM), 

In fact, a stronger result holds. 
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Proposition 1.1 The operators t and n can he extended to continuous surjec- 
tive maps 

t:H{d,n''M) H-^l^{d,Vl^dM), 
n: H{5,n^+^M) H-^^^{d,n^-''dM), 

where the space H^^/^{d,fl''dM) is the space of k -forms uj^ on dA-I satisfying 
This result is due to Paquet 

The formula (^6|) can be used also to define function spaces with vanishing 
boundary data. Indeed, let us define 

for all Ty'^'+i G H{S,n''+'^M)}, 
for all Lj'' e H{d,itM)}. 

It is not hard to see that indeed 

H{d, n'^M) = t-^{0}, H{6, = n-^{0}. 

We are now in the position prove the following lemma. 
Lemma 1.2 The adjoint of the operator 

d : L^in'^M) D H{d,^tM) ^ L^i^t+^M) 

is the operator 

and vice versa. Similarly, the adjoint of 
is the operator 

d : L^{n''M) D H{d,n^M) L^{n^+^M) 

Proof: We prove only the first claim, the other having a similar proof. 

Let rf''^^ G T>{d*), where d* denotes the adjoint of d. By definition, there exists 
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for all uj'^ e H{d, il'^M). In particular, if w*^ G fi'^M'"*, we see that, in the weak 
sense, 

i.e., 6r]''+^ = e L'^{VL^M). Thus, Ty'^+i e i7((5, O'^+^Af), and the claim follows 
now since we have 

for all w'^ G if(d,n'=M), i.e., (5 = d*. □ 

In the sequel, we will write for brevity H{d) = H{d, fl'^M), etc. when there is 
no risk of confusion concerning the order of the forms. 

For later reference, let us point out that the Stokes formula for the complete 
Maxwell system can be written compactly as 

(?/, Muj)l2 + {M-q, = (tw, n?7) + {t-q, nw) , (17) 

where w G H with 

H = H{d) X [H{d) n H{5)] X [H{d) n H{5)] x H{5) (18) 

and 77 G H^(Af) and we use the notations 

tu = {tuj^ ,tuj^ ,tuj'^) nuj = (nw^, ncj^, ncj-^), 

and, naturally, 

(tw, nrfj = {tuj°, nri^) + {tu^ ,mf) + (tcj^ nyy^). 

With these notations, we give the following definition of the Maxwell operators 
with electric and magnetic boundary conditions, respectively. 

Definition 1.3 The Maxwell operator with the electric boundary condition, de- 
noted by 

M,:ViM,)^L^{M), 

is defined through the differential expression (|J^, with the domain 2?(A^c) C 
L^(Ai") defined as 

V{M,) = H t := H{d) x [H{d) n H{5)] x [H{d) n H{5)] x H{5). 
Similarly, the Maxwell operator with the magnetic boundary condition, denoted 

by 

Mn. : V{Mrn) L2(M), 

is defined through the differential expression (|7^, with the domain I?(A^,„) C 
L^(M) defined as 

V{Mn,) = H „ := H{d) X [H{d) n H{S)] x [H{d) n H{S)] x (5). 
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Before further discussion, let us comment the boundary conditions in terms of 
physics. For vectorial representations of the electric and magnetic fields, the 
electric boundary condition is associated with electrically perfectly conducting 
boundaries, i.e., nx E = 0, n- B = 0, where n is the exterior normal vector at the 
boundary. In terms of differential forms, this means simply that t^;'' = tw^ = 
and t *o B'' — tw^ — 0. On the other hand, the magnetic boundary conditions 
represent a magnetically perfectly conducting boundaries, i.e., nx H = 0, n-D = 
0, which again in terms of forms reads as tH^ ~ t(l/a)i/^ = or t(l/Q;) * o;^ = 
nuj^ = and t *o = t{l/a)v^ = 0, or in terms of uj^ , t(l/a) * lu^ ~ nuj^ ~ 0. 

There is an obvious duality between these conditions. It is therefore sufficient 
to consider the operator with the electric boundary condition only. This obser- 
vation is related to the well-known Maxwell duality principle. 

Consider the intersections of spaces appearing in the domains of definition in 
the previous definition. Let us denote 

H\{n^M) ^ {uj^ e H^nHl) I tuj'' ^ 0}, 
h],{^^M) = {uj^ e H\n''M) I nuj'' = 0}. 
It is a direct consequence of Gaffney's inequality (see [|T]) that 

Hid, n'^M) n H{6, n'^M) = Hlin'^M), 
H{d, n'^M) n H{6, nHi) = niin^M). 



The following lemma is a direct consequence of Lemma 1.2 and classical results 
on Hodge- Weyl decomposition |Q. 

Lemma 1.4 The electric Maxwell operator has the following properties: 
i. The operator Aic is skew-adjoint. 

ii. The operator M.c defines an elliptic differential operator in M*"*. 

iii. Ker{Mc) = {{0,uj\uj'^,uj^) G H „ : duj^ = 0, Suj^ = 0, doj'^ = 0, Suj'^ = 
0, Slu^ = 0}. 

iv. Ran{Mc) = L^ifl^M) x {6H{6,n^M) + dHid,n°M)) x 
{SH{S, n^M) + dH{d, n^M)) X dH{d, Q^M). 

By the skew-adjointness, it is possible to define weak solutions to initial-boundary- 
value problems needed later. In the sequel we denote the forms uj{x,t) just by 
uj{t) when there is no danger of misunderstanding. 

Definition 1.5 By the weak solution to the initial boundary value problem 
Lut + Mlu = p e Ll^im,L^M)), 



12 



t'^laMxM = 0' '^(•,0)=a;oeL2, (19) 

we mean the form 

u{t) = U{t)ojo + [ U{t- s)p{s)ds, 
Jo 

where U{t) = exp(— tA^e) is the unitary operator generated by M.^- 

In the analogous manner, we define weak solutions with initial data given on 
i = T, T e K. Assuming p e C(K,L^(M)) and using the theory of unitary 
groups, we immediately obtain the regularity result 

uj e C(R, L2)nCi(M,H'). 

where H' denotes the dual of H. 

We shall need later the boundary traces of the weak solution. To define them, 
let (wow) Pn) S 'D{M.e) X C(M,'D(A4e)) be an approximating sequence of the pair 
{uJo,p) in \? X C(M,L2). We define 

U)n=ll{t)ijJon+ / U{t - s)pn{s)ds, 

Jo 

whence w„ e C{R,V{Me)) n C1(M,L2). Let if = {ip" , (p'^ , (p'^) be a test form, 
ip^ G CQ°{R,Q^dM). Let ry be a strong solution of the initial boundary value 
problem 

r]t + J^T] = 0, 
tr] = ^, r/(-,0) = 0. 

We have 

{r]{T),L0n{T))i,2 = / dt{r],u!n)dt 
Jo 

T 

{{Mri,ujn) + {ri,MuJn))dt, 
and, by applying Stokes theorem, we deduce 

(??(r),C^„(T))L2 =- / {ip^HLOn). 

Jo 

Hence, wc observe that, when going to the limit n oo, the above formula 
defines noj = lim„^oo nci-',,, G I>'(M, where 

V'{dM) = V'{n"dM) X V'{n^dM) x V'{Q^dM). 

We conclude this section with the following result. 
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Lemma 1.6 Assume that the initial data ujq is of the form ujq = (0, Wq, Wq, 0), 
where Sluq = 0, dujQ = and we have p = 0. Then the weak solution lo of 
Definition satisfies also Maxwell's equations i^ld^, i-e., uj^ — and 

LO^ = 0. 

Proof: As observed in Remark 1, to and, in particular, satisfies the wave 
equation 

in the distributional sense, along with the Dirichlet boundary condition tw" = 0. 
The initial data for is 

^0(0) = = 0, 

and 

Hence, we deduce that also — 0. 

Similarly, oj^ satisfies the wave equation with the initial data 

^3(0) = ^3 = 0, 

and 

As for the boundary condition, we observe that 

tduj^ = tLOf - tduj^ = dttLu'^ ~ dtLu'^ = 0, 

corresponding to the vanishing Neumann data for the function *uj^. Thus, also 
= 0. □ 

1.2 Initial— boundary value problem 

Our next goal is to consider the forward problem and the Cauchy data on the 
lateral boundary dM x R for solutions of Maxwell's equations. Assume that w 
is a solution of the complete system. The complete Cauchy data of this solution 
consists of 

{tuj{x,t),nuj{x,t)), {x,t)edMxR+. 

Assume now that uj corresponds to the solution of Maxwell's equations, i.e., we 
have = and lo'^ = 0. Consider the Maxwell- Faraday equation in (]lO|), 

+ duj^ 0. 

By taking the tangential trace, we find that tui^ = —dtu^, and further, 

tuj^{x,t)=uj^{0)- [ dtLu^{x,t')dt\ xedM. 
Jo 
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Similarly, by taking the normal trace of the Maxwell- Ampere equation in (|T^) , 



Lol - Slu^ = 0, 



we find that nuj} — dnui'^, so likewise, 

nuj\x,t) ^nuj^{0) + [ dnuj^{x,t')dt' , x e dM. (20) 

In the sequel we shall mainly consider the case w(0) = 0, when the lateral 
Cauchy data for the original problem of electrodynamics is simply 

tco = (0,/,- / df{t')dt'), (21) 



ncu - (0,5, / dg{t')dt') (22) 
Jo 

where / and g are functions of t with values in U^dM. 

The following theorem implies that solutions of Maxwell's equations are solu- 
tions of the complete Maxwell system and gives sufficient conditions for the 
converse result. 

Theorem 1.7 Ifuj{t) e C(M,Hi) n C\R,L'^) satisfies the equation 

ujt+MLj = 0, t>0 (23) 

with vanishing initial data uj{0) ~ 0, and uj'^{t) ~ 0, to^it) — 0, then the Cauchy 
data is of the form (^-(^. 

Conversely, if the lateral Cauchy data is of the form (pj|j-(p^ for < t < T , 
and LO satisfies the equation with vanishing initial data, then uj(t) is a 

solution to Maxwell's equations, i.e., w°(t) = 0, uj^{t) — 0. 

Proof: The first part of the theorem follows from the above considerations if we 
show that LL}{t) is sufficiently regular. 

Since u;^ e C(M, H^{n^M)) we see that nu;^ e C(M, H^^^{n^dM)) with dnu>^ e 
l,H-^/^{n^dM)). Furthermore, as Suj}{t) = dSuj^{t) = 0, 

ULul e C{R,H-^/^{n^dM)), Slu^ e C{]S.,H-^/^{n^dM)), 



which verifies (22). 

To prove ( |2l] ) we use Maxwell duality: Consider the forms 

a 

Then rj — (jy*^, t]^ ,r]^ , rf) satisfies Maxwell's equations rjt + Airj = where A4 is 
the Maxwell operator with metric g and scalar impedance a^^. In sequel, we 
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call Maxwell's equation with these parameters the adjoint Maxwell equations 
and the forms 77^ the adjoint solution. Now the formula (22) for adjoint solution 
implies (21) for uj. 

To prove the converse, it suffices to show that uj'^{t) = 0. Indeed, the claim 
uj^{t) = follows then by Maxwell duality described earlier. From the equations 

u'^-Slo'^ = 0, (24) 
ul +duj° -duj^, = (25) 

it follows that lu'^ satisfies the wave equation 

UJ^t + SdLU° = 0. 

It also satisfies the initial condition (jj"(0) = and (^t{0) = and, from (pTf), 
boundary condition toj° = 0. Thus, w° = for < t < T. 

□ 

The following definition fixes the solution of the forward problem considered in 
this work. 

Definition 1.8 Let f = if,f,f) e C°°([0, T]; f2(5M)) be a smooth bound- 
ary source of the form i.e., f^ — 0, ff — —df^. Further, let R be any right 
inverse of the mapping t. The solution of the initial-boundary value problem 

LUt+Mu) ^0, t > 0, 

w(0) = wo e L2(M), tuj^f, 

is given by 

Lu = Rf +U{t)LUo~ f U{t~ s){MRf{s) + Rfs{s))ds. 
Jo 

We remark that the boundary data / could be chosen from a wider class / S 
H^/^{dM X [0,r]). 

Theorem ll.Tl motivates the following definition. 



Definition 1.9 For solution uj of Maxwell equations we use the fol- 

lowing notations: 

i. The lateral Cauchy data for a solution uj of Maxwell's equations with van- 
ishing initial data in the interval < t < T is given by the pair 

{tuj^{x,t),nLj^{x,t)), {x,t) £ dM X [0,T]. 
ii. When uj satisfies initial condition uj{0) = the mapping 

: CS^{[0,T],n\dM)) ^ Co°§([o,r],r!i(aAf)), 

Z^(tw^) = nuj^\dMx[o,T], 
is well defined. We call this map the admittance map. 
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Here C^{[0,T], B) consists of C°° functions of t with values in a space B, i.e. 
B = ri^((9Af) in definition 1.9, which vanish near t = 0. 

Note that in the classical terminology for the electric and magnetic fields, Z'^ 
maps the tangential electric field n x E\qmxIo.t] to the tangential magnetic field 

n X H\gMx[0,T]- 

The boundary data and the energy of the field inside M are closely related. 
The following result, crucial from the point of view of boundary control, is a 
version of the Blagovestchenskii formula (see Q for the case of the scalar wave 
equation). Observe that the following theorem is formulated for any solutions of 
the complete system, not only for those that correspond to Maxwell's equations. 



Theorem 1.10 Let uj and 77 be smooth solutions of the complete system 
Then the knowledge of the lateral Cauchy data 

{tuj, nuj), {tr], nri), 0<t<2T, 

is sufficient for the determination of the inner products 

(L^^(t),V(s))L2, 0<i<3, 0<s,t<T 

over the manifold M . 



Proof: The proof is based on the observation that, having the lateral Cauchy 
data of a solution uj, we also have access to the forms dtu and dnw at the 
boundary. On the other hand, t commutes with d so that 

tdu!-' — dtuj^ , nSu!-' — t * *d— * uj^ — dt — * uj^ — dnoj^ . 

a a 

Let us define the function 

F^s,t) = {u;i{s),r^Ht)). 
From the complete system, it follows that 
id',~df)F^is,t) = iu;iM,ri'it))L^-iu^'is),4tit))L^ (26) 
= -{{dS + 6d)cj^{s),rjJ{t))L2 + iui^{s),{d6 + 6d)r]'it))L^ 
= fis,t). 

By applying Stokes theorem we obtain further that 

f^{s,t) = (nw^tV) + (tuj^,ndT]^) - (t^w^nV) - (nduj^ ,trf) , 

where we suppressed for brevity the dependence of the boundary values on s 
and t. Now the complete system implies that 
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and, similarly, 

df]-' ~ —iji^^ + Srf^'^, Syf = rjl ^ + dr/-'^'^. 

A substitution to the above formulas then gives 

P{s, t) = {nuj^,tr]t^ + dtT]^-^) + {tuj\ -nr]{+^ + dm]^+^) 

-(twf-i +c^W~^n^y^■) - (-nuji+^ + dnLj^+^ ,trf) , 

where d stands for the exterior derivative on DM. Hence, is completely 
determined by the lateral Cauchy data. What is more, we have 

(0, t) = (s, 0) = 0, F^^ (0, t) = Fi (s, 0) = 0. (27) 

Hence, we can solve F{s, t) using (p6| ) and ( [27| ) as claimed. □ 

Remark 3. If w and rj are solutions to Maxwell's equations (|o|-|ll|), the for- 
mulas above simplify. We have 

and 

Then, for j — I the inner product {uj^ (t) , uj^ (t)) 1^2 defines the energy of the 
electric field. Similarly, for j = 2 the inner product {u'^ {t) , oj'^ (t)) 1^2 defines the 
energy of the magnetic field. 

2 Inverse problem 

The main objective of this chapter is to prove the following uniqueness result 
for the inverse boundary value problem. 

Theorem 2.1 Given dM and the admittance map Zt, T > Sdiam (M), for 
Maxwell's equations, (fF^j- ([IJj, it is possible to uniquely reconstruct the Rie- 
mannian manifold, {M,g) and the scalar wave impedance, a. 

Observe that, once we know the travel time metric g as well as the wave 
impedance a, formula gives the metrics and g^, which correspond to 
the material parameters fi and e. 

The proof of the above result is divided in several parts. The first step, which 
is discussed in the next sections, is to prove necessary boundary controllability 
results. These results are used, in a similar fashion as in to reconstruct 

the manifold and the travel time metric. 
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2.1 Unique continuation results 



In the foUowing lemma, we consider extensions of_difFerential forms outside the 
manifold M. Let T C 9A/be open. Assume that M is an extension of M across 
r, i.e. M C M, r C int(M) and DM \ T C DM. Furthermore, we assume that 
the metric g and impedance a are extended smoothly into M as g, a. In this 
case, we say that the manifold wit h sc alar impedance (M, g, a) is an extension 
of{M,g,a) acrossT. (See Figure 2.1). 



Figure 1: Manifold M is obtained by gluing an "ear" to M. 
We have the following simple result. 

Lemma 2.2 Assume that M is an extension of M across an crpen set T C dM . 
Let u}^ he a k-form on AI and uj^ he its extension hy zero to M. Then 

1. Ifuj'' e H{d,n''M) and tcj'']^ = 0, then Zj^ e H{d,n''M). 

2. Ifuo^(^H{5,Vt^M) andnuo\^0, then U;'' e H{5,n''M) . 

Proof: External differential, in terms of distributions, of lj'^ can be defined by 

where ip'' S is arbitrary. However, by the formula (|l^), 

= (a;^(5¥.'=)i2(M) = (^^5^/)i2(M) + (ta;^n/). 

Moreover, since supp{(p'') CC Af'"*, then supp(n(p'^) C F, where tcu'' vanishes. 
Thus, 

i.e. duj'' is the zero extension of du''. In particular, duj^ E L'^{M), so that 
2'^ e H{d,^tM). 
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The claim concerning the codifFerential is proved by a similar argument. □ 
As a consequence of this result, we obtain the following. 

Theorem 2.3 Lef G C^K, L^) H C(M, H), tcj|rx[o,T] =0, ntj|rx[o,T] =0, be 
a solution of the equation ujt + M.oj = in M x [0, T] . Let lo he its extension 
by zero across T C dM . Then the extended form, w{t) satisfies the complete 
Maxwell's system on (M,g, S), i.e. uit + Mxo — in M x [0,T]. 



We are particularly interested in the solutions of Maxwell's equations. The 
following result is not directly needed but wc have included it, since the basic 
idea is useful when we will prove the main result of this section. 



Lemma 2.4 Assume that lu in the above theorem satisfies MaxwelVs equations, 
i.e., = and uj^ = 0, and w(a;,0) = 0. Iftu^ = and nw^ ^ on T x [0,T], 
then u) satisfies Maxwell's equations in the extended domain M x [0,r]. 



Proof: From Theorem 1.7 it follows that, since lo satisfies Maxwell's equations. 



tw = (0,tw\-/' dtuj'^dt') ^Q, nw=(0,nw^/" dnuj'^dt') ^ 
Jo Jo 

in r X [0,T]. Therefore, the previous theorem shows that the continuation by 
zero across F, uj{t), satisfies the complete system in M x [0,r]. 

However, w"(t) = 0, w^(t) = in M x [0,T], i.e., uj{t) satisfies Maxwell's 
equations with vanishing initial data in the extended manifold AI . □ 

When we deal with a general solution to Maxwell's equations, (p^-(pTl), which 
may not satisfy zero initial conditions, and try to extend them by zero across 



F, the arguments of Lemma 2.4 fail. Indeed, if ujq ^ 0, then ( pOD show that 
ncj^ = is not sufficient for nuj^ = 0. However, by differentiating with respect 
to time, the parasite term nti;^(0) vanishes. This is the motivation why, in the 
following theorem, we consider the time derivatives of the weak solutions. 

Denote by r(x, y) the geodesic distance between x and y on (M, g). Let F C dM 
be open and T > 0. We use the notation 

K{T, T) = {{x, t)eM X [0, 2T] \ t{x, T) < T - \T - t\} 



for the double cone of influence with base on the slice t = T. (see Figure 2.1 



Theorem 2.5 Let Lu{t) be a weak solution of Maxwell's system in the sense 
of Definition with cuo = (0,Wq,Wq,0). Assume, in addition, that 6ljq = 0, 
dujQ = and p ~ 0. Lf nuj^ ^ in Fx]0,2T[, then dtto = in the double cone 
K{T,T). 
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Figure 2: Double cone of influence. 



Proof: Let tp £ C^([— 1, 1]), J^^ ip{s)ds = 1 be a Friedrich's moUifier. Then, for 
any a > and uj{t) G C(]0, 2T[), L^(M)) satisfying conditions of the Theorem, 
denote by uj^{t) its time-regularization. 

Then € C°°{[a, 2T~a[, L^(Af)) continue to be weak solutions to the Maxwell 
system and, moreover, to Maxwell's equations (p^-(pl|). Thus, 

Mlu, = -dtoJa e C°°([a,2T- f7[,L2(M)), 

i.e. uja & C°°([ct, 2r — a-[,I?(Alo))- Repeating these arguments, 

As nuj„ = V'ct * nw, 

n{uj„f =0 on r X [a, 2T - (t[. 
Applying (pO[), we see that ndtWa = on F x [a,2T — a[. 

Denote by uj the extension by zero of cu across F and 770- that of dtUJa- We claim 
that, in the distributional sense, r]^ satisfies the complete Maxwell system, for 
a < t < 2T - a. Indeed, let (p = {ip°, (p^,(p^, (p^) € C^Qcr, 2T - a[, r2Af'"*) be 
a test form. Using the brackets [ • , • ] to denote the distribution duality, that 
extends the inner product 
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we have 



[dtija + Mri„,Lp] = -[i]„,Mf + (ft] 

As tuja = 0, it follows from the Stokes' theorem and the fact that uja satisfies 
Maxwell's equations, that 

j.2T I-2T 

Jo Jo 

As supp(t(y9) C rx](7, 2r — a[, where nujcr — 0, the right side of this equation 
equals to 0. In addition, tuia- = for t £]a,2T — (t[, where t is the tangential 
component on dM. Thus, the claim follows. 

However, rja G C°° (]cr, 2T — cr[, i^))- Therefore, similar considerations to the 
above shows that this implies that 

GC°°([a, 2T-a[,V'^{M,)), 

i.e. rja is infinitely smooth in Af"^* x [a, 2T — a[. Since rj„ = outside M x E, the 
unique continuation result of EUer-Isakov-Nakamura-Tataru , that is based 



on result of Tataru 1 55 



for smooth solutions, implies that i]a- = in the double cone t{x, M \ M) < 
T — a — \T — t\, X G M, where t is the distance on (M, g). As rjcr = dtWa in 
M, this implies that dtuj^ = in the double cone 

r(x,r) < T-CT- |T-<|, xeM. (28) 



When a —I- 0, rja- ^ dtto, in the distributional sense, while the cone (|2§| ) tend to 
K{r^ T) and the claim of the theorem follows. □ 

We note that the unique continuation result of is related to scalar e, /i. 
However, it is easily generalized to the scalar impedance case due to the single 
velocity of the wave propagation. 

Following the proof of Theorem 2.5, we can show the following variant of The- 
orem 
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Corollary 2.6 Let Lo{t ) be a weak solution to the complete Maxwell system in 
the sense of definition |i.4 with p — 0, and, in addition, on rx]0,T[. // 

T > 2 diam{M), then uj'^{t) = 0, '^'^{t) = and uj{t) is a solution of Maxwell's 
system for < t < T . 

Proof: We will consider only cj" using the n Maxwell duality for u^. By remark 
1 and (H), 

c^o + = 0, tc^°|9Mx[o,T] = 0. (29) 
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Also 

+ dLO° - 5uj^ = 0, 

imply, together with (p2[), that 

xiduj^ = nSuP' — nuj] = rfncj^ — nuj] = 

on r X [0, T]. Together with the boundary condition in (p9|), this shows that 
the lateral Cauchy data of uj{t) vanishes on F x [0,r]. Using now the wave 
equation in (p9|), this imply that, due to Tataru's unique continuation |]55| , 
Wo = in the double cone K{r,T). As T > 2diam(M), this yield that 
wO(T/2) = c^?(r/2) = 0. It now follows from (||) that Lu°{t) = for < i < T. 
□ 



2.2 Introduction for controllability 

In this section we derive the controllability results for the Maxwell system. We 
divide these results in local results, i.e., controllability of the solutions at short 
times and in global results, where the time of control is long enough so that the 
controlled electromagnetic waves fill the whole man ifold. Both types of results 
are based on the unique continuation of Theorem 2.5 and representation of inner 
products of electromagnetic fields over M, in a time slice, in terms of integrals 
of the lateral Cauchy data over the boundary dM over a time interval which is 



given by Theorem 1.10 



Consider the initial boundary value problem 

ujt+MLj = 0, t>0. 



(30) 



with the initial data a;(0) — and the electric boundary data of Maxwell type, 



tio = (0, /, 



dfit')dt') 



(31) 



where we assume that / e C^{R+,ft^dM). By Theorem 1.7, we know that 
uj°{t) = and uj^{t) = 0. 



Let uj denote the weak solution of Definition 1.5 with p — and w(r) — loq. 
Assume, in addition, that the conditions of Lemma 1.6 are satisfied so that uj 
satisfies also Maxwell's equations. 

As we have seen, Stokes formula implies the identity 



{lj{T),luo) = ~ f {tuj,nuj)dt. 
Jo 



(32) 



We refer to this identity as the control identity in the sequel. 



23 



2.3 Local controllability 

In this section, we study differential 1— forms in M that can be generated by 

using appropriate boundary sources active for short periods of time. Instead 
of a complete characterization of these forms, we show that there is a large 
enough subspace in L^(f2^M) that can be produced by boundary sources. The 
difficulty that prevents a complete characterization is related to the topology of 
the domain of influence, which can be very complicated. 

Let r c dM be an open subset of the boundary and T > arbitrary. We define 
the domain of influence as 

M(r, T) = {xgM\ t{x, r) < T}, 

where r is the distance with respect to the travel time metric g. Observe that 

M(r, T) = ^^(r, T)n{t = t}. 

Furthermore, let oj be the strong solution of the initial-boundary value problem 

wt+Muj = 0, a;(0) = 0, 

with the boundary value 

ta; = (0,/,- tdf{t')dt'), 
Jo 

where / e Cg^QO, T[, Qip) with C^{]0,T[,n^T) being a subspace of forms in 
Cq°(]0, T[, Q^dM) with support in T. To emphasize the dependence of co on /, 
we write occasionally 

(j = ujf = (0,(a;■^)^(w■^)^0). 

We denote 

X{T,T) = clL.{{cof)'{T) I / e C^{]0,T[,Q'T)}, (33) 

i.e., X{T, T) is the L^-closure of the set of the electric fields that are generated 
by CQ°-boundary sources on rx]0,T[. Furthermore, we use the notation 

H{5, M(r, T)) = {w^ e H{5, M), supp (a;^) e M(r,T)}. 

We will prove the following result. 

Theorem 2.7 The set X{T,T) satisfies 

6H^{n'^M{r, T)) c X{T, T) c c1l2 [5H{5, M{T, T))^ . 

Here Hl{^l'^S), 5 C M is a subspace of HQ{Vt^M) of forms with support in 
cl(5). 
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Proof: The right inclusion is straightforward: Since to satisfies Maxwell's equa- 
tions, we have = and, hence. 



uj\T)^ [ Suj^dt £ 6H{d,M{T,T)). 
Jo 



To prove the left inclusion, we show that any field of the form = Si]'^ with 
7?^ e H^{M{T,T)) is in (X(r,r)-L)^. To this end, let us first assume that 
Wq e X^(ri^M) is a 1-form such that 

for all uj^ = i^^HT))'^ generated by boundary sources / G {]0,T[,n^r). 
Since ui^ = 5liP , it suffices to consider only those forms Wq that are of the form 
Wq = (5?7q for some ?7o G -fflt^)- Indeed, by Hodge decomposition (see [^) in 
L2(M), we have 

where dSp = 0, tSp = so Wq _L tj^ automatically. 

Let w be a weak solution, at the time interval [0, T], of the initial boundary 
value problem (|l^) with tuj — 0, and 

5(.,r) = (o,c^o\o,o) = wo. 

By our assumption, 

and thus, by the control identity, ( ^ ) and conditions 

Jo Jo Jo 

for all differential 1-forms / G Co°°(]0, r[, f^^r). Thus, we have 

nS^ = on rx]0,T[. 

Furthermore, it is easy to see that, for T + t G [T, 2T], we have 

^{T + t) = {0,w\T - t),-^\T - t), 0), 

and, therefore, also 

nw^ = on rx]T,2T[. 

But this implies that, as a distribution, nw^ vanishes on the whole interval 
]0, 2T[ since it is in ^^^i^, H-^/'^{dM)). By applying the Theorem [2^, we can 
deduce that if = in the double cone K{T,T). In particular, we have that 
dujl = Zjj(T) = in M(T,T). 
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Let now = 5rf € SH^{n'^M{T,T)). Then 

{i'^,ujI)l2 = {rj^,dh)l)L^ = 0. 

This holds for arbitrary wl S X{T, T)^, i.e., v e (X(r, T)^)^ = X{T, T). □ 

Remark 4. Later in this work, we are mainly interested in controlling the time 
derivatives of electromagnetic fields. Let us denote 

C°°(r,T) = { f f{t')dt' I / G C^{%T[,n^T)}. 
Jo 

With this notation, we have 

XiT,T) = cli.{(a;/(T))i | / G C~(r,T)}. 
Indeed, if = {oj^^ G X{r,T), then {ujf)^ = (wf)!, where 



m= f f{t')dt'. 
Jo 



J, u 

Conversely, the time derivative of a field cj-' , / G C°°(r, T) satisfies the initial- 
boundary value problem with the boundary source ft G Cq"(]0, T[, fi^F). 

2.4 Global controllability 

We start by introducing some notations. Let u) be the strong solution of the 
initial-boundary value problem 

U!t+MLU = 0, w(0)=0, 

with the boundary value 



ta; = (0,/,- / df{t')dt'), 
Jo 



where / G Cg°{]0, To[, O^L), Tq > and T c dM is an open subset. For T > Tq, 
we define 

Y{T, T) = {LoliT) I / G C^QO, Tol fi^F)}. (34) 

For F = dM we denote Y{T) ~ Y{dM,T). Our objective is to give a charac- 
terization of the set Y{T) for Tq large enough. In the following, 

rad(M) = maxT(a;, dM). (35) 
We prove the following result. 
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Theorem 2.8 Assume that To > 2rad(M). Then, for T > Tq, clL2(Af)F(T) is 
independent of T , i.e. cli,2(^pj^Y{T) =Y, and, moreover, 

r = {0} X SH{S) X dHid) X {0}. (36) 

Remark 5. The result holds also for Y replaced with Y(r,T), when 

Tq > 2maxT(a;,r). 

Proof: Let w = w-^ be a solution. As / = for T > Tq, we have 

tuj\T) = 0, 

and, consequently, 

u;t{T) = -Mlj{T) = {0,5Lj'^{T),^dLj\T),0) 
e {0} X dH{d) X dHid) X {0}. 

To prove the converse inclusion, we show that the space Y{T) is dense in {0} x 

6H{6) X dHid) X {0}. To this end, let wq e {0} x x dHid) x {0} and 

luq _L YiT). This means that, for arbitrary lo = lo^ satisfying the initial- 
boundary value problem ([l^), 

(cc.o,^t(T))L2 = (c^i,c^i(r))i2 + iu;l,uliT))L2 = 0. (37) 

Let Zj denote the weak solution of the problem 

u)t + Muj = 0, 

tuj = 0, w(r) — Wo- 

Observe that the initial value satisfies 

SujI = 0, dujl = 0, 

which implies that uj satisfies Maxwell's equations. Consider the function F : 
K ^ R, 

Fit) - (ii(i),Wt(i))L2. 

We have, by using Maxwell's equations, that 

Ftit) = (w,tJ4t)L2 + (5i, Wi)L2 

— —(5"'", iduj^)i^2 — (o;^, dSuj^)]^2 + iduj^, duj^)i^2 + iSuj'^, Suj'^)]^2, 
and further, by using Stokes' theorem, 

Ftit) = -^{tu^it),ndLo\t)) - (n52(t),t5cj2(t)). 
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However, to; = and SuP' = lo}. Thus, 



On the other hand, the initial condition lo{{)) = 0, together with the orthogo- 
nahty condition @, imply that F(0) = F{T) = 0, so that 

/ {n^^Jt)dt^-~ [ Ft{t)dt = 0. 
Jo Jo 

Since / £ Cg°{]0,Tin^T) is arbitrary, this implies that 

nl^jf = in rx]0,r[. 



But now Theorem 2.5 implies that cutt vanishes in the double cone K{T,T/2). 
By the assumption Tq > 2rad(M), this double cone contains a cylinder 

C = Mx]T/2~ s,T/2 + s[ 



with some s > (See Figure 2.4). Therefore, lou that satisfies Maxwell's 




Figure 3: The double cone contains a slice {T/2} x M and the waves vanish 
near the slice t = T/2. 

equations and homogeneous boundary condition tujtt — 0. Therefore, it vanishes 
in the whole M x M. In particular, this means that, with some time-independent 
forms wi and ^2, 
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with t(jji = 0, tcj2 = 0. Again, by Maxwell's equations, we have 
for all t. Therefore, 

UJ2 = MlUi, MuJ2 = 0. 

But then. Stokes' theorem implies that 

(w2,W2)l2 = {UJ2, MuJl)i,2 = -{MuJ2,1^i)l2 = 0, 

i.e., 0J2 — and Aiuji — 0. Observe that, by the assumption of the Theorem, 
= S(r) = Wo = (0, -Si^^, di^^,0) = Mv, 

o 

for some v G {0} x H{d) x H{S) x {0}. Therefore, a further application of Stokes 
theorem gives 

(tJi, wi)l2 = {uJi,M:^)-L2 = -{Muji, v)-L2 = 0, 
i.e., also ui — loq = 0. The proof is therefore complete. □ 

2.5 Generalized sources 

So far, we have treated only smooth boundary sources and the corresponding 
fields. For later use, we need more general boundary sources. 

Let Y = c1l2(m)^(9M, r) be the space of the time derivatives of electromag- 
netic fields satisfying Maxwell's equations, see (p^, (^). We define the wave 
operator 

: Co°°(]o, r[, n^dM) ^y, / ^ c^/(r), 

where T > Tq and Tq > 2rad(M). By means of the wave operator, we define 
the J-"— norm on the space of boundary sources as 

11/11^ ^IIW^'^/IIl- (38) 

The definition of this norm is independent of the choice of T > Tq by conserva- 
tion of energy. 



Notice that by Theorem 1.10 , the knowledge of the admittance map Z^-^ enables 
us to calculate explicitly the J^-norm of any smooth boundary source. 

To complete the space of boundary sources, let us define the equivalence ~ of 
sources by setting 

/ ^ g iff W^f = W^g. 
Further, we define the space J^([0,To]) as 

.F([0,To])=C7o°°(]0,r[,r!i5M)/~. 
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Finally, we complete T{[0,Tq]) with respect to the norm (J3^). Hence, this 
space, denoted by JF([0,To]) consists of Cauchy sequences with respect to the 
norm (p8[), denoted as 

/-(/,)}to, f,eC^{]0,T[,n'dM). 

Note that, for any f G we can find h G !F such that h — f and h = {hj)j^i, 
hj_G C^(]e,T[,rjiaM) for some e > 0. The reason for this is that Theorem 



2.8 is valid also with Tq replaced with Tq — e, when e is small enough. Thus, 
for small £ > 0, we can define, for any f = (/j)°^o ^ ^^^'^ translation 

/(•+^) = a(-+£)),°^oG^- 

These sources are called generalized sources in the sequel. The corresponding 
electromagnetic waves are denoted as 

uj{{t) = lim u;l'{t) for t > Tq. (39) 

By the isometry of the wave operator, the above limit exists in for all gen- 
eralized sources. 

We note that the above construction of the space of generalized sources in well- 
known in PDE-control, e.g. ||9|, H]. 

Remark 6. Observe that since the wave operator is an isometry and 
J([0,To]) was defined by closing C^(]0, T[, rj^aM) with respect to the norm 
(p8|), the wave operator extends to a one-to-one isometry 

7^ Lo{iT), T{[o,To]) ^ ciL2(r(aM,r)), 

where the target space is completely characterized in the previous section. 
We say that h E T is a generalized time derivative of / e JF, if for T ^ Tq, 

Jirn^\\t±^l^^H\\^= (40) 



= lim II ^ - (T)IIl.(m) = 

In this case we denote h — B/, or just h — dt f- In the following, we use 
spaces ^F" = 'D{W), s £ Z+, which are spaces of generalized sources that have 
s generalized derivatives. Note that, if (^) is valid for T = Tq, it is valid for 
all T > Tq due to the conservation of L^-norm for Maxwell's equation (energy 
conservation). Thus, if / G .F*, we have, for T >Tq, 

M^u({T) = dl4{t)\t=T e L^M). (41) 



Note that A4 here is the differential expression given by (|15|), rather than an 

4 



operator with some boundary conditions. Since tuj{ ^ on dMx]To, oo[, we 
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see that t(a/wf)(t) = for t > To and j <s- 1. Thus, for / e J^" and T > Tq, 

_ S 

4 ^ f]{C'~^ {[T,oolV{Mi)) nR&n{Me)). (42) 



Moreover, by ([4 21) and Lemma 1.4 



a;/(r) G Hf„,(Af™*) for T > Tq. 

Next we consider dual spaces to the domains of powers oi Aie- Since C 
I?(7W^), we have (I?(A^c))' C H-^ Similarly, we see that Hg " c (2?(7W^))'. 
These facts will be needed later in the construction of focusing sources. 

2.6 Reconstruction of the manifold 

In this section we will show how to determine the manifold, M and the travel 
time metric, g from the boundary measurements of the admittance map Z. 
We will show that the boundary data determines the set of boundary distance 
functions. The basic idea is to use a slicing principle, when we control the 
supports of the waves generated by boundary sources. 

We start by fixing certain notations. Let times Tq < Ti < T2 satisfy 
To > 2 rad(M), Ti > To + diam(M), T2 > 2 Ti 

We assume in this section that the admittance map Z"^^ is known. 

Let Tj C dM be open disjoint sets, I < j < J and t~ and Tj' be positive times 
with 

< r" < r+ < diam(Af), I < j < J. 
We define the set S = S{{Tj,T^,T^}) C M as an intersection of slices, 

,7 

S=f] {M{T„T+)\M{r„T-)). (43) 
i=i 

Our first goal is to find out, by boundary measurements, whether the set S 
contains an open ball or not. To this end, we give the following definition. 

Definition 2.9 The set Z = Z({rj, t^^, r^^})^^]^ consists of those generalized 
sources f € J-°° {[0,Tq]) that produce waves ujt = U}( with 

1. LuliT,) e X{T^,T+), for all j, 1<J<J, 

3. uJttiTi) = in M(r,, Y), for all j, l<j<J. 
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M 



Figure 4: The sources / of the waves uj^{x, t) are supported on the time- interval 
[0, To] which enables us to control the waves at times t > Tq. In the construction 
of the manifold, the supports of the waves arc considered at time t = Ti. To 
this end, we use the unique continuation in double-cones (triangle in the figure) 
intersecting the boundary in the layer dM x [To,Ti]. Note that in this layer it 
is crucial that / = 0. 

Remark 7. Observe that, since uit satisfies Maxwell's equations, we have, in 
particular, 

These identities imply that, at t = Ti, ivtt is of the form 

LOttiTi) = {0,0,ojUTi),0) = {0,0,d'n\0), 
for the 1-form r]^ = — , and 

supp((i77^) C S. 

This observation is crucial later when we will discuss focusing waves. 
The central tool for reconstruction the manifold is the following theorem. 

Theorem 2.10 Let S and Z be defined as above. The following alternative 
holds: 

1. If S contains a,n open hall, then dim(Z) = oo, 

2. If S does not hold an open ball, then Z = {0}. 

In order to prove the above alternative, we need the following observability 
result that will be also useful later. 

Theorem 2.11 Given the boundary map Z^'^ , we can determine whether a 
given boundary source f e .F°°([0,To]) is in the set Z or not. 
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Figure 5: In Definition 2.£ we can consider e.g. the case Fi = T, = si, 
Ti = 0, and F2 — dM, T2 — diam (Af), — S2- Then the waves that 

satisfy the definition have the following properties: By 1., the wave {uj{)^{Ti) 
coincides with a wave that is supported in A/(F, si). This domain of influence 

on the figure is the upper part of the cone of influence. Thus, d{uj{Y{Ti) = 

(a;/()^(ri) is supported in M(F, si). By 2., the wave (w/j)^(ri) vanish in the 

boundary layer M{dM, 52)- Combining these, we see that lo(^{Ti) is supported 
in A = M(F, si) \ M{dM, S2). 



Proof: Let / = {fk)kLo •^°°([0,7o]) be a generalized source. Consider first 
the question whether (lljI)^{Ti 
to the existence of a sequence, 



the question whether (lljI)^{Ti) G X(Fj,r^+). By Remark 4, this is equivalent 



such that 



JmJ{cjt)\T,)-{u;'^^)\r;)\\=0. 



(44) 



By linearity of the initial-boundary value problem, we have 

\\iu;t)\T^)-iu;'^^)\r+)\\^\\iu;<'--n\T^)l 

where the source gk.e is 

9kAt) = ilk)t{t) - {hMt + r+-T^) e c^{%T^[-n^dM). 

However, by Lemma l.lCl| , || (a;^'= *)^(Ti)|| is completely determined by the ad- 
mittance map, Z'^'^ making possible to verify 



In a similar fashion. Condition 2 of the definition of Z is valid for /, if 



hm ||(^/'=(Ti))^|| =0, 
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and this condition can also be verified via the admittance map, Z,'^'^ . 

Finally, consider Condition 3. We assume here that we already know that / 

satisfies Conditions 1-2. First, we observe that iott = (^{t satisfies 

{dt+M)uJtt = in M x M+, 

along with the boundary condition 

tutt in 9M X [To,oo[. 

If Condition 3 holds, by the finite propagation speed, Utt vanishes in a double 
cone around Tj, i.e., 

uJtt = in Kj = {{x, t) e M X M+ | r(x, Fj) + |i - Ti] < t"}, 

for all J = 1, . . • , J ■ In particular, this means that, in each Kj, ut does not 
depend on time, and Condition 2 implies that ujf = in Kj. Hence, we have 

nuj^ = Z'^V on F^- x ]Ti - r^" , Ti + [. (45) 

Conversely, assume that condition (|4|) holds together with Conditions 1-2. 
Then ujt satisfies 

{dt+M)uJt = in Af X M+ 
with the boundary conditions 

tLol = 0, nujf = in F^- x ]Ti - t" , Ti + r" [. 

Here we used the fact that Ti — tJ' > Tq^so that / = in Fj x ] Ti — tJ' , Ti + t^" [. 
Now the Unique Continuation Principle, given by Theorem |2.5| , implies that 
LOtt = in Kj and, in particular. Condition 3 is valid. The proof is complete as 
it is clear that the condition (^) is readily observable if the admittance map, 
2^T2 , T2 > Ti + Tj , is known. □ 

Now we can give the proof of Theorem 2.1C. 



Proof of Theorem 2. It: Assume that there is an open ball B d S. Let 
7^ (y5 € be an arbitrary smoot h 2- form with supp ((/?) C B. From the 

global controllability result. Theorem 2.8, it follows the existence of a general- 
ized source / G J^{[0, To]) such that 

u;/(ri) = (0,(5^,0,0). (46) 
Moreover, ip £ il^B implies that ip £ P(A^g) for any s > so that / £ 

^°°([o,ro]). 

We will now show that f £ Z. Indeed, Conditions 1-2 are obvious from the 
definition (^) of /. Finally, we observe that 

4(Ti) = -MujfiTi) = (0,0,-d^^,0), 
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so Condition 3 is also satisfied. This proves the first statement of the theorem. 

To prove the second part, assume that S does not contain an open baU. Suppose, 
on the contrary to the claim, that there is a non-vanishing source f € Z which 
produces the wave uj{t) — uo^ {t). Then, by Conditions 1 and 2 in Definition 2.9 



Furthermore, 
so that 



J 

supp(u;t(ri)) c fl M(r„T+) = s+. 



On the other hand, Condition 3 in Definition 2.9 imply that 

cottiTi) = in U/=i M{T„Tp = . 

Thus suppa;t((Ti) C \ . However, if the set S does not contain an open 
ball, then the set 5*+ \ 5*^ in nowhere dense. Since LOtt{Ti) is smooth in Af'"*, 
it vanishes in M . In particular. Maxwell's equations imply that 

duj]{T{)^-ujl{T,)^0. (47) 



On the other hand, ujt e c1l2 (y(5M, Ti)), so Theorem |2j implies that lo]{Ti) 
is of the form 

for some 2-form r/. Since, in addition, tLol{Ti) = it then follows, by Stokes 
formula, that 

{loI{T,),uj]{T^))l2 = = (77^da;l(^l))i2 = 0. 

Together with Condition 2, this implies 

contradicting to the assumption / 7^ 0. The proof is complete. □ 

We are now ready to construct the set of the boundary distance functions. For 
each X S M, the corresponding boundary distance function, is a continuous 
function on dM given by 

: dM ^ R+, r^iz) t{x, z), z e dM. 

They define the boundary distance map TZ : M C{dM), TZ{x) — r^, which is 
continuous and injective (see Q , |^ ) . We shall denote the set of all boundary 
distance functions, i.e., the image of TZ, by 

7^(^f) ^ {r^ e C{dM) I X e M}. 
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It can be shown (see Q, [£3) that, given TZ{M) C L°°{dM) we can endow it, 
in a natural way, with a differentiable structure and a metric tensor Tj, so that 
{TZ{M),g) becomes an isometric copy of {M,g), 

{n{M),g)^{M,g). 

Hence, in order to reconstruct the manifold (or more precisely, the isometry 
type of the manifold), it suffices to determine the set, TZ{M), of the boundary 
distance functions. The following result is therefore crucial. 

Theorem 2.12 Let the admittance map Z'^^ be given. Then, for any h G 
C{dM), we can find out whether h G TZ{M). 

Proof: The proof is based on a discrete approximation process. First, we observe 
that the condition h S Tl{M) is equivalent to the condition that for any sampling 
Zi, . . . ,zj € dM of the boundary, there must be x S M such that 

h{zj) — t{x, Zj), ^ < j < J- 

Let us denote Tj — h{zj). By the continuity of the distance function, t{x,z) 
in a; G M, z S 9M, we deduce that the above condition is equivalent to the 
following one: 

For any e > 0, the points Zj have neighborhoods Vj C dM with diam(rj) < e, 
such that 

int ( fl M(r, , r, + e) \ M{T, , r, - e)] ^ 0. (48) 



On the other hand, by Theorem 2.1C, condition (ttSh is equivalent to 



dimi Z{{Tj, Tj + £, Tj — e} I —00 



that, by means of Theorem 2.11, can be verified via boundary data. □ 



As a consequence, we obtain the main result of this section. 

Corollary 2.13 The knowledge of the admittance Z^^ is sufficient for the re- 
construction of the manifold, M endowed with the travel time metric, g. 

Having the manifold reconstructed, the rest of this article is devoted to the 
reconstruction of the wave impedance, a. 

2.7 Focusing sources 

In the previous section it was shown that, using boundary data, one can control 
supports of the 2-forms {uj(^f'{t). In this section, the goal is to construct a 
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sequence of sources, (/p), p — 1,2, ■ ■ ■ such that, when p ^ co, the corresponding 

forms (wf^)^(Ti) become supported at a sing le point, while {uj([y{Ti) = 0. For 
t>Ti, these fields behave like point sources, a fact that turns out to be useful 
for reconstructing the wave impedance. 

In the following, let S denote the Dirac delta at y g M'"*, i.e.. 



Sy{x)<p{x)dVg{x) = <p{y), 

M 



where g C^{M) and dVg is volume form of {M,g). 

Since the Riemannian manifold {M, g) is already found, we can choose Tjp C 

3P 3P 



DM, < T p < T+ < diam(Af), so that 



dSp, f]Sp^{y},ye M^"*. (49) 



0=1 



Then, Zp — Z{{Tjp, Tj ,T f}'^^^}) is the corresponding set of generahzed sources 



defined in Definition 2.9 



Definition 2.14 Let Sp, p = 1, 2, • • • , satisfy (]44^. We call the sequence f = 
(/p); P — ■ • • with fp G Zp, a focusing sequence of generalized sources of 
order s (for brevity, focusing sources), s G Z_|-, if there is a distribution-form 



A = Ay on M such that 

lim {io^*^-{Ti),fj)^2 = {Ay,fj)^2, 

p—^oc 

for allrie V{Ml). 

Remark 8. Observe that, by the identity, 

Luf- = 4" (50) 

and Remark 7, the electromagnetic wave cjf'"'^^ (Ti) is supported in cl(S'p), so Ay 
must be supported on {y}. 

We will show the following result. 

Lemma 2.15 Let the admittance map Z'^'^ be given. Then, for any s S 2+ and 
any sequence of generalized sources, (fp), p = 1,2 ■ ■ ■, one can determine if (fp) 
is a focusing sequence or not. 

Proof: Let r/ £ T>{M^). We decompose as ry = r^i + 772 , where 

r]ieV{Ml)nc\{Y), ?72 e2?(Mc)nr^. 



37 



By the global controllability result, Theorem 2.S, and isometry of the wave 
operator W^,T>To, 

771=4, /Je^^([0,To]). 

Since Lof*^'' € clL2(y), so that lo^^'^'' -L r]2, the condition that / is a focusing 
source is tantamount to the existence of the limit 

iAy,r,)= lim(t^f^''(ri),4)i., (51) 

p — *oo 

for all h e jr*([0,To]). However, by Theorem 1.1C| , the existence of this limit 
can be verified if we are given Z'^^ . 

Conversely, assume that the limit (|l]) does exist for aU h e J^''([0, To]). Then, 
by the Principle of Uniform Boundedness, the mappings 

?7K^ (a;f'-^-(Ti), 77)^2, p e Z+, 

form a uniformly bounded family in the dual of T>{Ai^). By the Banach-Alaoglu 
theorem, we find a weak*-convergent subsequence 

cjf'-{T,)-.Aye(v{Ml)] , 



which is the sought after limit distribution-form. □ 

Since supp(Aj,) is a point, Ay consists of the Dirac delta and its derivatives. The 
role of the smoothness index, s is just to select the order of this distribution, as 
is seen in the following result. 

Lemma 2.16 Let Ay = limp^^ ujf'^" {Ti) is a distribution of order s = 3. 
Then Ay is of the form 

Ay{x) = iO,0,d{XSy{x)),0), (52) 

where X is a 1-form at y, AG T*M. Furthermore, for any A € T*M there is a 
focusing source f — with Ay of form ^5^. 

Proof: From the results in Section |2.5| , we deduce that, when s = 3, 

Ay e {V{Ml)y c H-3. 
Furthermore, from Remark 7, the electromagnetic waves ( |50| ) are of the form 

c^f^''(Ti) = (0,0,rfr,p,0), (53) 



for some 1-forms rjp. Combining ( these with the fact that supp(^y) = {?/}, we 
^,0). Here Ay, expressed, for exampl 
) near y, must be of the form 

Al{x)^a^S(x)e,+V''dkSJx)e,, 



see that Ay = (0, 0, A^, 0). Here Ay, expressed, for example, in Riemann normal 
coordinates {x^,x'^,x ) near y, must be of the form 
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where 6j = {l/2)ejki.dx^ f\dx^ and is the totally antisymmetric permutation 
symbol. Furthermore, by (p3|), 

dAl = {a^dj5y{x) + V''dkd,5y{x))dVg = 0. 

Let if he a, compactly supported test function and, in the vicinity U of y, 

(p{x)=x\ 1,2,3. 

It follows that 

Further, let -0 be a compactly supported test function and, in the vicinity U of 

y, 

i;{x)=x^x'', j,fc = l,2,3. 

As before, we obtain 

0= {dAl,^j) = b''' + b''^. 
Thus, V'' may be represented as b^'^ — e^'^^Xi, \i G T*M , implying that 

Al{x) = e^''Xedk5y{x)e,. (54) 

By the properties of the permutation symbols, e^*^^, 

e^'=^6lj = ieJfc^e^.p^rfa;P a dx« = S^S^dxP A 



Substituting this expression back to (|54|), we finally obtain 

Alix) = \idk5y{x)dx^ A dx'^ = d{6y{x)\idx^), 

as claimed. □ 

By the above results, for any y S Af™* and A G T*M , we can, in principle, find 

focusing sequences / such that uj(^{Ti) = Ay, where Ay is of form (p2|). We 
should, however, stress that, at this stage, we can not control the corresponding 
\ = \{y). 

Consider now a family of focusing sources fy, y d M'"* , with the corresponding 
1— forms X(jj). 

Lemma 2.17 Given the admittance map Z^^ , it is possible to determine whether 
the map y i— > Ay is a nowhere vanishing 1-form valued C°° -function. 

Proof: Let ip e fi^M'"* be an arbitrary compactly supported test 1— form. By 
Theorem 2.8, there is a generalized source h S J-""" such that 
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Let / — {fp),p = 1, 2, • • • , be a focusing source of order s = 3. Then, by Lemma 
2.16 and the definition of the focusing sources, we have 

lim(4(Ti),c.~^ri)) = {Ay,J^iT,)) (55) 

p — >oo 



JM 

Further, by Maxwell's equations, 

AyA*v'(y), 

Here, for A, t? € T*M, 

XA*ri= (A, ■r])y dVg = g^^Xj-qk dVg. 

Thus, 

hm {4-{T^),u^{T,)) = (A, ^{y))y. (56) 



By Theorem LIO , the inner products on the left side of equation (|56| ) are obtain- 
able from the boundary data. Thus, we can find the map y {X,ip{y))y, y €E 
M'"*. Since ip € rj^M™* is arbitrary, this determines whether A e fi^M'"*. It 
also determines whether Ay = or not for any y G M™*. This yields the claim. 
□ 



Another way to look at Lemma 2.17 is that the admittance map, Z^'^ determines, 
for any boundary source h G C^{]0,To[;^^dM), the values, at any y € M'"*, 
of (A, {u;^t)^{t))y for some unknown A G ^I^M'"* and Ti < t < T2 - diam(M). 
Moreover, using this map, we can verify that the 1— forms Xk{y), corresponding 
to three families of focusing sources fk{y), k — 1,2, 3, are linearly independent 
at any y g M'"*. These give rise to the following result. 

Lemma 2.18 Let Z'^'^ be the admittance map of the Riemannian manifold with 
impedance {M,g,a). Then, for Ti < t < T2 — diam(M) and h G !F{]0,To[), it 
is possible to find the forms 

L{y)U\y,t))\ K{y){u^iy,t)f, 

at any y € A/'"'. Here L{y) : T*M ^ T*M and K{y) : A^T*M A^T*M are 
smooth sections of End{T* M™^) anc? End(A^r*M'"*) , correspondingly. 



We emphasize that, at this stage, L{y) and K{y) arc unknown. However, they 
are independent of t or h. 

Proof: As M is already found, we can choose three differential 1— forms, G 
rj^M"^* which, at any y G Af"^*, form a basis in T*M. Using the families fkiy) 
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of focusing sources introduced earlier, we can construct, for any h E !F°°{]0, To[), 
the differential 1— form, 

p^{y,t) := (Afc(y), 6(y). (57) 

This defines a smooth section, L{y) of End(T*Af'"*), 

L{y){u;}{yM ^ P~\y,t) en'M''^', 

proving the assertion for (a;^(y,i))^ with h € J-"°°(]0, ro[). Its extension to 
h e ^(]0, To[) is an immediate corollary of the fact that ^°°(]0, ro[) is dense in 
J^{]0,To[) in the J^-norm. 

To analyse {ujf{y,t))'^ , consider the form 

7] ^ (— *iLj^, — *uj'^, — *uj^, — *uj°) (58) 
a a a a 

(cf. J/- forms in formulae (^) and This form satisfies the complete Maxwell 

system ^ 

rjt + Mrj ~ 0, 

where M. is the differential expression (p^), corresponding to the manifold 
{M,g,a~^). Then the admittance map t?7^|aMx]o,T2[ ^ ^a-^v'^\dMy.]n.T^[ is 
the inverse of the given admittance map Z'^'^ : tw^|gMx]o.T2[ ~^ ^a^'^\dM x]o,T2[- 
Thus, Z^'^ determine the admittance map, for M and we can apply the 

results for {lo^Y to {rjlY , where / = Z^^h. Namely, we can find L{y){rit{y,t)y , 
where L{y) : T*M T*M is a smooth section of End(T*Af'"*) which, at this 
stage, is unknown. At last, since 

*L{y){rj{{y,t)y=K{y){u;}{y,t))\ 

for some smooth section, K{y) of End(A^r*M'"*), the assertion follows. □ 

We note for the further reference, that, similar to the case of the 1— forms, the 
construction of K{y) involves a choice of three differential 2— forms, that we 
denote by /Xfe € fi^Af™*, and three families of generalized sources Kkiu) that 
satisfy 

wf'^^(ri) = (0,%fe^^(a;)),0,0), f,keA^T;M, fc = 1,2,3. (59) 

Below we call the generalized sources fk (y) the focusing sources for 2-forms and 
Kfc(j/) the focusing sources for 1-forms. 

Before going to a detailed discussion in the next sections of the reconstruction 
of a, let us explain briefly the main outline of this construction. It follows from 
Lemma |2.18 that, using the admittance map Z, we can find the electromagnetic 
waves ujj (t) , Ti < t < T2 — diam (Af), up to unknown linear transformations, L 
and K. We observe that, by Theorem for any basis ^fe(j/), A: = 1, 2, 3, there 
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arc families fkijj) of focusing sources, such that the corresponding transforma- 
tion L is just identity. Indeed, to achieve this goal, we should choose fk{y) in 
such a manner that, at any y G M™*, 

icjty^riT,) = diXk{y)Q, (60) 

where Xk{y) is dual to C/c(y), 

{My),^ji.y))v=^kj- (61) 

In the next sections we will identify conditions on fk{y) and Kk{y), verifiable in 
terms of Z, which make L and K to be identities. 

2.8 Reconstruction of the wave impedance 

In the previous section, it was shown how to select a family fy, y £ M'"* of fo- 
cusing sequences such that the corresponding electromagnetic fields concentrate 
at f = Ti at a single point, y, 

lim ui-{Ti) = {0,0,d{XSy),0). 

Here X = Xy G T*M is yet unknown. Moreover, we can select this family in 

such a manner that, as a function of y, Xy G Jl-'^M™'. In particular, it means 
that for times < t < ty = T{y, dM), the electromagnetic wave defined as 

Ge{y) = Ge{x,y,t) = G4X]{x,y,t) = lim u{i{t + Ti), (62) 

satisfies the initial-boundary value problem 

{dt+M)G^{y) = OmMx]0,tyl 

tGe(y) = OmdMx]0,ty[, (63) 

Ge(2/)|t=o = (0,0,d(A4),0). 

The solution to this problem is called the electric Green's function. We will use 
this solution to reconstruct the scalar wave impedance, a on M. We start with 
analysis of some properties of Ge- To this end, we will represent Gg in terms 
of the standard Green's function, G = G{x,y,t) for the wave equation on 1— 
forms. Thus, G is defined as the solution to the following initial-boundary value 
problem 

(5^ + dS + Sd)G{x, y, t) = {d^ + Al)G{x, y, t) = X5y{x)5{t) in M x R 

tG{x,y,t)= 0, (64) 
G{x,y,t)\t<o = 0, 
where A e T*M is a given 1-form. 

This Green's function has the following asymptotic behaviour. 
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Lemma 2.19 For < t < ty, Green's function G{x,y,t) for the 1-form wave 
equation (jfe^, has the representation 

G{x, y, t) = 6{t - t{x, y))Q{x, y)X + r{x, y, t). 

Here Q{x,y) : T*M — > T*M is a bijective map that corresponds to a (1,1)- 
tensor depending smoothly on {x, y) G Af'"* x M'"* \ diag{M™'^). The remainder 
r{x,y,t) is a bounded function, when t < ty, where ty is small enough. 

The proof of this lemma is postponed in the Appendix. 

In the following, we fix y e Af'"' and X = Xy E TyM . By operating with the 
exterior derivative d on the both sides of the differential equation in (|6j), we 
see that 

{d^,+^l)dG{x,y,t) = d{X5y)5_{t). 

Hence, using the decomposition {d'j + A) = [dt + M){dt ~ M.) , we find that the 
form ui{t) — (0, 0, dG{x, y, t), 0) satisfies the equation 



(dt + M) i^idt - M)iO, 0, dG{x, y, t),0) j = Dy^xm, 

where 

Dy^x = iO,0,d{Xdy),0). 
Let Gc{y) — Gc{x,y,t) be defined as 

G,{x,y,t) ^ {dt-M)iO,0,dG{x,y,t),0) 

= iO,SdGix,y,t),dtdGix,y,t),0). 

Then, due to the finite propagation speed, G{y)\Qj^^^^^ ^ j = 0, so that Gc{x, y, t) 

satisfies the boundary condition tGdy) — for t < ty. Invoking the uniqueness 
of solution for (|63|), we see that G^iy) = Gdy), t < ty. 
Now, using Lemma [2.19 , we obtain that 

dG\{x,y,t) = {Q{x,y)Xy A dT{x,y))5}^\t - T{x,y)) +ri{x,y,t), 

where S}-^"^ is derivative of the delta-distribution and the residual ri is sum of 
the delta-distribution on dBy (t) , where By (t) is the ball of radius t centered in 
y, and a bounded function. Thus, we see that 

dtdGl{x,y,t) ^ {Q{x,y)Xy/\dT{x,y))5j'^\t-T{x,y))+r2{x,y,t), 
SdGl{x,y,t) = *{dT{x,y)A*{Qix,y)XyAdT{x,y)))S^^\t-T{x,y)) 

+r3{x,y,t), (65) 

where residuals r2 and r^ are sums of first and zeroth derivatives of the delta- 
distribution on dBy(t) and a bounded function. Moreover, by formulae (|65|), 

^By(t)[Q{x,y)Xy AdT{x,y)]^Q, (66) 
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nB„(t) [* {dT{x, y) A *(Q(x, y)\y A dT{x, y)))] = 0, (67) 

where tBy{t)^^ , T^By{t)^^ are the tangential and normal components of lo^ on 
dByit). This corresponds to the physical fact that the wavefronts of the electric 
and magnetic fields are perpendicular to the propagation direction. Let now 
fy be a focusing source for a point y e M'"*. Due to the definition (|T]) and 
the definition of the generahzed source (p^), there is a sequence {f^)^^i, € 
C^(]0,To[;r2iaAf), p = 1,2,- •• such that 

J^{T^) = hm uihTi), 

p — *oo 

and the right side is understood in the sense of the distribution- forms on JIM'"*. 
Then, for t > 0, 

wf" (i + Ti ) = lim Lu(^ (t + Ti). (68) 

p — >C30 

Applying Lemma ^.18 , it is possible to find, via given Z^^, the magnetic com- 
ponents K{x){uj{ {x,t + Ti))^, f — fjj oi these fields with K being a smooth 
section of End(A2r*M'"'). At last, using (||), we find 

K{x){Lu{''{x,t + Ti)f = lim K{x){JJ {x,t + Ti)f e V {n^ M'"^). 

p— i-OC 



Since 

we see that 

when t <ty. In particular, we can find the singularities of Green's function up 
to a linear transformation K{x). 

Hence we have shown: 

Lemma 2.20 Let fy = {f^), p = 1,2,- ••, be a focusing source for a point 
y. Then, given the admittance map, Z'^'^ , it is possible to find the distribution 
2— form KGc{y)^ for all x satisfying T{x,y) < ty, where ty is small enough. In 
particular, the leading singularity of this form determine the 2— form 

K{x){Q{x,y)XyAdT{x,y)). (69) 



As shown at the end of the previous section, K E End(A^r*M"^*) was obtained 
by using three focusing sources Kk{x), k = 1,2, 3. Our next goal is to formulate 
conditions, verifiable using boundary data, for K to be isotropic, i.e. 

K{x) = c{x)I, (70) 
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Figure 6: Vector v is the right singularity of the electromagnetic wave in the 
plane M x {t}. The reconstructed singularity w has wrong direction, if the 
transformation matrix K{x) is not isotropic. 



with c{x) being a smooth scalar function. 

We start with observation that, for a given Kk{x), k — 1,2, 3, and any fy, we can 
find tB^{t)K{-){Q{-, y)AAdr(-,?/))|^, for x e dBy{t) and smaU t > 0. Here K{x) 
is the linear transformation corresponding to the chosen Kk{x). This follows 



from Lemma 2.20 and the fact that the underlying Riemannian manifold (M, g) 



is already found. When K is isotropic, it follows from (|66|) that 

tB^it)K{■){Q{■,v)\^dT{■,y))=0. (71) 
Let us show that the condition (|7l|), which is verifiable via boundary data. 



actually guarantees that K is of form (|70|). 

Indeed, for a given A G T*M and any x, t with sufficiently small t — t{x, y), ( [n]) 
means that K{x){Q{x, y)X/\dT{x, y)) continues to be normal to the 2— dimensional 
subspace, T^{dBy{t)) C T^M, i.e. 

[K{x){Q{x,y)XAdT{x,y))]{X,Y) = 0, X,Y & {T,dBy{t)). 

As Q{x,y) : TyM —> T*M is bijective, Q{x,y)X A dT{x,y) runs through the 
whole 2— dimensional subspace in A^T*M, normal to Tx{dBy{t)), when A runs 
through TyM. Therefore, (|7l|) implies that K{x) keeps this subspace invariant. 

Let us now vary y and t keeping x fixed. Then Tx{dBy{t)) runs through the 
whole Grassmannian manifold, 6*3^2 (TxM). Therefore, ( |7l| ) implies that, if 
S A^T*M is normal to a subspace C € Gz^2^ then K{x)liP remains to be 
normal to C. This is implies that the eigenspace of K{x) has dimension three 
and hence we must have 

K{x) = c{x)I, 

where c(x) is a scalar function and / is the identity map in A?T*M . 
In the following, we always take focusing sources which satisfy (|70|). 
Thus, we can find the values of the 2-forms 

{^l{x,T^f =def c{x)(J,{x,T^)f, / e T. 
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Our further considerations are based on the equation, 

d{4{Tif) = 0. (72) 

Thus, we intend to choose those focusing sources Kk{x) which produce K{x) = 
c{x)I with 

In this case, 

= d{c{x){4T{)f) = dc{x) A (^/(Ti))2, 

due to By the global controUabihty, the form {lj{)'^{x, Ti) runs through the 
whole A^*M. This imphes that dc{x) = at any x S Af'"'. Hence, c{x) — Cq 
is a constant. Thus, we choose focusing sources, Kk{x), so that 

{^}{T^)f = coiu;}{T,)r. 
Evaluating the inner products, 

/ wUx,Ti)A*i:;fix,Ti)^ [ c^cjUx,Ti) A *coUjUx,Ti), 
Jm Jm 

we can compare them with the energy integrals, 

-^ujUx,T,)A*uj^{x,T,), 



which can be found from the boundary data by means of Theorem 1.10 . By 
considering waves ui^ {x,Ti) with 

lim supp(w/'(-,ri)) = {y}, 

wc find the ratio 

/^#(a:,ri))^A*(^f-(x,ri))^ , 
lim — ~ %a{y). 

The above considerations imply that, using Z'^^ , we can determine a{x) up to 
a constant Cg. Since the impedance map satisfies ZM,g,ca = c~^ZM,g,ca, we see 
that, knowing the impedance map, we can also determine cq. 

As the fact that Z'^ , T > 4 dia mAf determines {M,g) is already proven, this 
completes the proof of Theorem 2.1. □ 
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2.9 Back to W 



In this section we use the obtained uniqueness result for Maxwell equation on a 3- 
diniensional manifold to analyze to the group of transformations which preserve 
the boundary data in the dynamical inverse problem for Maxwell's equation 
in a domain C M''. 

Assume that two Maxwell systems with electric and magnetic permittivities 
e'^(x), nl.{x) and e^(cc), /i|,(a;), x G il, have the same admittance map on 
dVl X [0, r], where T is sufficiently large. Denote by {M,g,a) and {M,g,a) the 
corresponding abstract Riemannian manifolds with impedance. By Theorem 



(M,5,a) = iM,g,a) 

i.e., there is an isometry H : {M,g) — > {M,g) and a — H^,a.. We can represent 
the abstract manifold [M, g, a) as the domain with the metric tensor, gij{x), 
given in Euclidean coordinates by (|^), and scalar impedance, a(x) in these 
coordinates. Similarly, we represent manifold {M,'g,a) using Euclidean coordi- 
nates in and obtain the metric tensor Ijij and impedance a. Then H : M —t M 
corresponds to a diffeomorphism, 



X-.n^n, X\9n = tdlan, (73) 



and 



g = X^g, i.e., g^'{x) = ——gP'^ix), x ^ X{x), (74) 

a — X^,a, i.e., a{x) — a{x). 

Using ( [t^ and (^), we see also that 

— f)x^ f)x-^ 

g, = X,g,, i.e., gl^{x) = ——gP'^{x), x = X{x). (75) 

Employing formula (|^), we obtain 

^l = ^egr^q, ^, = ^/f.gr^q- (76) 
Combining formulae (|7^)-(|76|) and introducing 



we obtain 



Similarly, 



Det {DX) dxP dxi 



^^.P9(:^), ^^X{x). (77) 
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Clearly, if X : fl ^ fl and X|go = id, the admittance map , T > is 



preserved in transformations (73), (|77|), (|78|) 



Thus we have proven the following result. 

Theorem 2.21 The group of transformations for Maxwell's equations ([^-(Q^ 
with a scalar wave impedance, which preserves the admittance map 2,"^ , T > 
4diam(M, 5), is generated by the group of difjeomorphisms offl satisfying (7c ). 
The corresponding transformations of e and fi are then defined by formulae (77), 



Remark 9. It follows from ([77|), ([78|), that e^'' and /x-''^ do not transform like 
tensors. This is due to the special role played by the underlying Euclidean metric 
= S^^ , which does not change by diffeomorphisms X. It should be noted that 
this form of transformations is observed also in the study of the Calderon inverse 
conductivity problem. Indeed, it is shown in that, for the conductivity 
equation in C M^, the boundary measurements determine the anisotropic 



conductivity up to same group of transformations as described in Theorem 2.21 
The Calderon problem is closely related to the inverse problem for Maxwell's 
equation, for instance, the low-frequency limit of the admittance map Z°° is 
related to the Dirichlet-to-Ncumann map for the conductivity equation |38| . 



2.10 Outlook 

There are several direction to which the present work can be extended. 

1. Natural question is the minimal observation time required to parameter 
reconstruction. It can be shown that the admittance map for any t > 0. 
Thus, it follows from the above, that Z"^ , T > 2radAf determines uniquely 
the manifold M, metric g and wave impedance a. The reconstruction of Z* 
for any t > may be obtained by a direct continuation of the admittance 
map, i.e., without solving the inverse problem. This continuation is a direct 
generalization to the considered Maxwell's case of the technique developed in 
[ pT| , p5|] for the scalar wave equation. An analogous method has recently be 
applied to Maxwell's equations in |p^ . 

2. Another natural inverse boundary value problem is the inverse boundary 
spectral problem for the electric Maxwell operator Me defined in Definition 



1.3, The problem is to determine the metric g and wave impedance a, or, 
in the other words, e and /x from the non-zero eigenvalues Xj of Me and the 
normal boundary values of the corresponding eigenforms. This problem was 
studied in, e.g. [|6|, for scalar Maxwell's equations. For the considered 
anisotropic case, this requires significant modifications of the method developed 
in this paper and will be published elsewhere. 

3. It often occurs in applications that the measurements are made only on a 
part of boundary. In formalism of this paper this means that we actually know 
only the restriction of the admittance map to the part F x [0, T] of the lateral 
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boundary, i.e. we are given 

For the scalar wave equation the corresponding problem is studied in Q , Q 
(see also [^). The combination of these methods and those of 

the present paper will be useful for analyzing the corresponding problem for 
Maxwell's equations. 

Appendix: The WKB approximation 

Here we consider asymptotic results for Green's function or, more precisely, for 
Green's 1-form, G{x,y,t) — G\{x,y,t), which is defined as the solution for the 
wave equation 

d^Gx + {dd + 6d)Gx=aSy{x)d(t) mMxR+, (79) 
Gx{-,y,t) = fort<0, tGx{-,y,t) = 0, 

where 6 ^ 6a. Here, A is a 1-form A = ^^^iKdx^ in normal coordinates 
{By{p),X), X = (x^, x^, x^), near a point y € M'"*, X{y) — 0. We assume that 
By{p) n dM — 0. In these coordinates, Sy{x) — S{x), when x £ U. Clearly, we 
can find, instead of the solution to (|7^), the fundamental solution 

d^G + {d6 + Sd)G = IS{x)6it) in M xR+, (80) 
G{-,y,t)^0 fort<0, tG(-, i) = 0, 

where / is the 3x3 identity matrix. Equation (|80|), written in normal coordi- 
nates, becomes a hyperbolic system 

{{d^ - g'^d.dj)! + B'd, + C}G = IS{x)S{t). (81) 

Here g^^ {x) is the metric tensor in these coordinates with 

g'^ (0) ^ S'^ , dkg'HO) ^ 0, (82) 

and B^{x), C{x) are smooth 3x3 matrices. We note that, in normal coordinates, 
T{x,y) = |a;|. However, we prefer to keep the notation t to stress the invariant 
nature of considerations below. 

Following ||l^, 1^, which deal with the scalar case, we search for the solution to 
(|l]) in the WKB form: 

G{x,t) = Goix)Sif~T^)+Y,Giix)Si-,it^-T^), (83) 

where Si{s) — s'_^/T(l + 1). Substitution of expression ( |83| ) into equation ( ^l| ) 
gives rise to the recurrent system of (transport) equations. The principal one is 
the equation for Go, 

dG 

4r_Jl(rJ) + {{g'Hrx) ^^^JT^{Tx) -6)1 + B\tx) d,T'^{Tx)) Go{tx) = 0, (84) 
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where x = x/t. In addition, to satisfy initial conditions IS{x)S_{t), corresponding 
to the right side in the wave equation (pO[), we require that 

Go(0) = ^I. (85) 

By (|8|), — 6 is a smooth function near a; = and g^^ didjT^\^_^ — 6 — 0. 

Also, dir'^l „ = 0. Therefore, 

' la:— ' 

^ { {g'^ (tx) didjT^ (tx) -6)1 + B' {tx) d^r^ (tx) } 

is a smooth function of (r, x) , so that Gq {x) is a smooth 3x3 matrix of x for 
r > 0. 

Actually, the matrix Go{x) is smooth everywhere in the neighborhood of y, 
i.e. X = 0, including a; = itself. Indeed, if we write the Taylor expansion 
of {g'^\x)didjT^ {x) — 6) / + B'^{x)diT^{x) near a; = and divide the result by 
T = \x\, we obtain that 

ig'^ {x)d^djT\x) ~ 6)I + B\x)d,T^{x) = ^ T^^rl^'l-ix''. (86) 

I/3|>1 

Substituting the Taylor expansion (with respect to r) of Go(t, x), 

Go(T,^) = ^Go,p(^)T^ 

into ( p4[ ) and using (|8^), (p5|), we obtain that Go,p(x)r^' are homogeneous poly- 
nomials of X of degree p: 

Go,p(J)tP = ^ Go,;3a;''. 

I/3|=P 

Then, Gg(x) = E |/3|<p ^o,/3 a: satisfies 

+ { (9^' 9,d,T^ - 6) / + 9,r2 } Gg = 0^, (87) 

where 

9P{t,x) = 6lP(a;) G G°°(C/), 6iP(a;) = ©(r^). (88) 
We construct Go as Gq (/ + Gq). Substituting this expression into (^) and 



using (BTp, (88), we obtain that 



where 



4T^(r, x) = AP{r, x) + yl^'(r, 3?) Gg(T, 5?), GP(0) = 0, 



AP{t,x) = -{GI{t,x))-^ 9P{t,x). 
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Therefore, Ap{t,x) S C°° as a function of {t,x) and Ap{t) = 0{tP). This 
imphes that Gq{t, x) = 0{t'p) and is C°° smooth as a function of (r, x), so that 
Gq, considered as a function of x = is in CP{By{p)). As p > is arbitrary 
and the solution Go of (||), (||) is unique. Go G C^{U). 

For G/, ^ > 1, we obtain transport equations 

4t^ + {(4^ - 6 + 9''{x)d,djT^ix)) I + B\x)d,T^} Gi 

and G/(0) = 0. If we write G; = GqFi, we obtain for Fi the equations 

+ 4ZFj = Go 1 [<?'^9,aj/ - B'd, - C] Gi^i, Fi[0) = 0, (89) 
with their solutions 

Fi{x) = f Go\sx) { [g''d,djl - B'd, - C] Gi^i] {sx) s'-^ds, (90) 

4 Jo 

being a smooth function of x. 

As (^) is a hyperbolic system, it is easy to show that the right side of ( p3[ ) 
represents the asymptotics with respect to smoothness of the Green's 1— form 
G{x,y,t), when t < T{y,dM). Clearly, ( ^3|) can also be written in the form 

G{x, t) = Go{x) 6{t^ - T^) + r{x, t) 

where r{x,t) is a bounded 3x3 matrix. 
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